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Abstract 

i-Ch ■ Several types of variable stars are found along the HR diagram whose pulsations 
are driven by the K-mechanism. Given their nature, the precise (T e s — L) 
O ' domain where these pulsators are located is highly dependent on the value 
of opacity and on its variation inside the star. We analyze the sensitivity of 



^ ■ opacity driven pulsators of spectral-type A and B (S Scuti, (3 Cephei and SPB 
<\ stars) to the opacity tables (OP/OPAL) and to the chemical composition of 
the stellar matter. We also briefly discuss the effect of opacity on pulsators 
whose oscillations are not driven by the K-mechanism, such as 7 Doradus and 
solar-like stars. 



Session: STARS - opacity driving, levitation, opacity data 
Individual Objects: 44Tau 

Introduction 



Already in 1919 Eddington proposed that an increase of the absorption coeffi- 
cient in some parts of the star at the moment of maximum compression could 
excite stellar pulsations. However, it was necessary to wait around forty years 
■ before Zhevakin (1959, and references therein) identified the region of the sec- 
\ ond ionization of He as responsible of pulsation excitation in cepheid variables. 
Baker & Kippenhahn (1962) made the first numerical stability computations 
applied to realistic models of stellar envelope and confirmed that the increase 
of opacity in the region of second ionization of He leads to sufficient excitation 
to overcome the damping in the rest of the star, explaining in this way the 
pulsation of Cepheids. Since then, other pulsators have been observed whose 
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excitation mechanism is linked to the variation of opacity with temperature and 
pressure, not only in the region of second ionization of helium, but also in the 
Fe-group opacity bump at logT ~ 5.3 {[3 Cepheids and SPB's eg. Dziem- 
bowski et al. 1993; sdB's, Kilkenny et al. 1997, Charpinet et al. 1996) and at 
logT ~ 6.25 (Wolf-Rayet, Townsend & MacDonald 2006). 

The necessary conditions for the excitation of a pulsation mode of period 
r by the K-mechanism were described by Pamyatnykh (1999): i) the opacity 
perturbation in the high temperature phase must grow outward; ii) the ampli- 
tude of the pressure eigenfunction must be large and change slowly with radius 
in the driving region; Hi) the thermal relaxation time (Airr 3 pC p T '/ 'Lr) in the 
driving region must be of the order or longer than the oscillation period (r). 
The behavior of stellar opacity as a function of temperature and density is not 
only important in the framework of the first condition. Stellar opacity deter- 
mines the density distribution of the equilibrium model and hence the oscillation 
frequencies. Moreover, the details of the stellar structure are in turn important 
in the balance between excitation and damping of oscillations. Finally, opacity 
of stellar plasma influences the location of convective regions, and hence may 
play a relevant role in the properties of oscillations. 

The main ingredients in the computation of stellar opacities are two: atomic 
physics, and chemical composition of the stellar plasma. In this paper we 
examine the consequences of recent updates on the oscillation properties of 
different types of main-sequence pulsators. 

Effect of opacity on stellar structure and evolution 

Opacity tables 

The improvements included in the new OP (Opacity Project) opacity compu- 
tations as well as a thorough comparison with OPAL (Iglesias & Rogers, 1996) 
opacity tables are described in Badnell et al. (2005). These changes result in an 
enhancement of Rosseland mean opacity, kr, at high density and temperatures, 
and an increase of 18% of opacity in the Z-bump due to the new Fe atomic 
data. The new OP opacity tables are much closer to the OPAL ones than the 
previous release (Seaton et al. 1994). Some differences remain, however, at 
low temperatures (logT < 5.5): the OP Z-bump in kr presents a hot wing 
slightly larger than the OPAL one. This translates into an kr difference that 
can reach 30% at low densities. On the other hand, OP and OPAL agree at 
high density within 5-10%. 

This disagreement should affect differently models of typical solar-like, 5- 
Scuti, and j3 Cephei stars. In figure 1 we plot OP and OPAL evolutionary tracks 
for models with typical values of mass for these three kind of variables (1.0, 
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Figure 1: Upper panels: evolutionary tracks of 1 Mq (left), 1.85 Mq (center) and 
12 Mq (right) computed with OPAL (solid lines) and OP (dashed lines) opacity 
tables, same chemical composition (X = 0.70, Z = 0.02) and treatment of convective 
transport. Lower panels: differences between OP and OPAL kr through an OPAL 
model at the middle of the MS (mass fraction of hydrogen at the center X c =0.35). 
Vertical lines indicate the boundary of convective regions in those models. 



1.85, and 12 M Q ). The 1.0 and 1.85 M Q OP evolutionary tracks are cooler than 
the OPAL ones. Since all the stellar parameters in the models are the same, 
these results suggest a larger OP opacity with respect to the OPAL one. The 
differences across the stellar structure of main sequence models are shown in the 
lower panels of figure 1. We emphasize a difference of ~ 30% at logT ~ 5.3 
for the 12 M Q model, that will have important consequences on the instability 
strip of B-type pulsators, and also that OP opacities are systematically 10% 
lower than the OPAL ones in the region around logT ~ 6. 



Solar mixture 

The re-analysis of solar spectrum by Asplund et al. (2005) (hereafter AGS05) 
including NLTE effect as well as tri-dimensional model atmosphere computa- 
tions, has led to a significant decrease of C, N, O and Ne solar abundances 
and as consequence, to a solar metallicity ((Z/X) & ) 30% smaller than the 



4 



Input from opacity data in computation of pulsation instability 



value provided by the "standard" GN93 mixture (Grevesse & Noels 1993). The 
abundance of iron (log N (Fe) / N (H)) has not been affected by the new anal- 
ysis, which means an iron mass fraction in AGS05 mixture 25% larger in than 
in GN93 one for a given Z. It is worth noting that while metallicity is mainly 
determined by the abundance of C, N and 0, the value of the Rosseland mean 
opacity can be strongly affected by other less abundant elements. The contri- 
bution to opacity from each element in Z depends on density and temperature, 
hence the effects of solar mixture changes on stellar pulsation properties are 
expected to depend on spectral type. Figure 2 shows the variation of kr with 
temperature through the structure of MS stellar models with masses of 1.0, 
1.85 and 12.0 M Q , and for each of them, the contribution of C, N, 0, Ne 
and Fe-group elements (Fe, Ni, Co, Mn) to kr. For 1.85 M Q (typical S Scuti 
star), kr (thin solid line) clearly shows the peaks due to H, and He ionization 
(logT ~ 4.1), to He+ ionization (logT ~ 4.6), and to the so called Z-bump 
of opacity at logT ~ 5.3). The opacity peak due to He + ionization is also 
evident in kr curve for 12.0 M© model, nevertheless, for this model the main 
contribution to opacity comes from Fe-group elements. We also notice that the 
contribution of CNO to kr decreases as stellar mass increases. 




log T log T lag T 



Figure 2: Profile of Rosseland mean opacity (thin solid line) for the same models as in 
figure 1. For each model we derived, by using the OP-sever facility 2 , the contribution 
to kr in percent coming from O (thick solid lines), N (dotted lines), C (short-dashed 
lines), Ne (dash-dotted lines) and from the elements of the iron group (long-dashed 
lines) 

The total effect of new solar mixture on the value of kr is shown in figure 3. 
It should be noticed that for a given value of metal mass fraction Z, AGS05 
mixture leads to a larger kr than GN93 one. Nevertheless, it should be kept in 
mind that the new Zq is only 70% of the old one. The hollow at logT ~ 6.2 in 
<5kr for 1.0 and 1.85 M (left and central panels in Fig. 3) is the signature of the 
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lower oxygen-abundance in AGS05 mixture. While the 30% decrease of C,IM,0 
abundances increases the discrepancies between the standard solar model and 
helioseismology, AGS05 values are in a better agreement with spectroscopic 
measurements in B-type stars (Morel et al. 2006). Moreover, since Fe-group 
elements are the main contributors to opacity for these stars, the increase of 
the iron mass fraction by 25% in the AGS05 mixture will favorably affect the 
excitation of (3 Cep and SPB pulsation modes in early-type stars. 




Figure 3: Relative differences of kr values provided by GN93 and AGS05 metal 
mixtures at fixed Z = 0.02, through the structure of the GN93 models. The models 
considered in left, central and right panel are the same as in figures 1 and 2. Vertical 
lines corresponds to the boundaries of convective regions. 



Opacity driven pulsators 

B-type pulsators 

(3 Cep and Slow Pulsating B (SPB) stars (spectral type B0.5-B3 the former, 
and B3-B8 the latter) are MS pulsators excited by the K-mechanism due to the 
Fe-group opacity bump at T ~ 2 x 10 5 K (e.g. Dziembowski et al., 1993). OP 
opacity values in the Z-bump are larger than OPAL ones by almost 30% for a 
typical (3 Cep, moreover the increase by 25% of iron mass fraction for a given 
Z implies a higher and larger peak in at T <~ 2 x 10 5 K. The combined 
effect of new opacity and solar mixture on the SPB and (3 Cephei instability 
strip has been examined by Miglio et al. (2007a, b) and Pamyatnykh (2007). 
Their results can be summarized as follows: while the different profile of n in 
OP and OPAL computations modifies the blue border of the instability strip, a 
larger Fe-mass fraction in the metal mixture provides slightly wider instability 
bands, and this effect increases as the metallicity decreases. Thus, the number 
of (3 Cep pulsators expected with AGS05 is more than three times larger than 
with GN93, and the SPB instability strip resulting from OP calculations is 
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3000 K larger than with OPAL ones. As a consequence, the number of expected 
hybrid j3 Cep-SPB objects is also larger for OP models. The Fe-mass fraction 
enhancement in the AGS05 mixture, compared with GN93, has the main effect 
of extending towards higher overtones the range of excited frequencies. While 
for Z = 0.01 the increase of B-type pulsators resulting from the updates of 
physics in the models is remarkable, for Z = 0.005 SPB-type modes are excited 
when considering OP with AGS05, but none of the different OP/OPAL and 
GN93/AGS05 evolutionary tracks for masses up to 18 M Q predicts (3 Cep 
pulsators. 

The updates in basic physics lead hence to decrease the disagreement be- 
tween theory and observations. In fact, the presence of B-type pulsators in 
low-metallicity environments (Kolaczkowski et al 2006, Karoff et al. 2008, Di- 
ago et al. 2008), as well as that of hybrid (3 Cep-SPB pulsators (Handler 
et al. 2006; Jerzykiewicz et al 2005) cannot be explained in the context of 
OPAL opacities and GN93 mixture. Nevertheless, some difficulties remain, for 
instance, on how to explain the presence of (3 Cep pulsators observed in the 
Small Magellanic Cloud (Kolaczkowski et al 2006), and the pulsation spectrum 
detected in v Eri (Dziembowski & Pamyatnykh 2008). 

The understanding of pulsations in subluminous B stars (sdB) has also 
benefited from the updated opacity values (Jeffery & Saio 2006), showing that 
B-type pulsators can be considered as a "critical test for stellar opacity" . 

S Scuti stars: 44 Tau 

5 Scuti stars (A-F spectral type) are located in the low luminosity end of the clas- 
sical instability strip. They pulsate in low order p and mixed g-p modes excited 
by the K-mechanism in the second ionization region of helium, at logT ~ 4.6. 
In a recent paper Lenz et al. (2008) have analyzed the effects of solar mixture 
and opacity tables in the seismic modeling of 44 Tau, a slowly rotating <5 Set 
variable that shows radial and non-radial pulsation modes (Antoci et al. 2007; 
Zima et al. 2007). The values of periods of fundamental (IIo) and first overtone 
(III) radial modes together with the low rotational velocity (v sin? < 5 km s _1 ) 
allow Lenz et al. (2008) to use the Petersen diagram (IIi/IIo vs log IIo) to de- 
rive the fundamental parameters of 44 Tau and to analyze their dependence 
on the basic physics. Comparison between OP and OPAL opacities through a 
typical 5 Scuti structure (central panel of fig. 1) does not show a significant dif- 
ference (lower than 5%) in the driving region. Nevertheless, the authors found 
that for the case of 44 Tau the OPAL opacities are preferable to the OP ones, 
since it was not possible to fit the radial modes of 44 Tau at the same time 
than its location in the HR diagram by using OP models. 

The Petersen diagrams (Petersen & Jorgensen, 1972) are known to be very 



J. Montalban, and A. Miglio 



7 




Figure 4: Relative differences between 
OPAL and OP kr trough the structure 
of an OPAL model with M = 1.875 M 
Z = 0.02, X = 0.70 and central mass 
fraction of hydrogen X c — 0.4. Dotted- 
line represents the relative difference of 
sound speed at fixed mass for the two 
models with the same parameters but dif- 
ferent opacity tables (OP et OPAL). 
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sensitive to the metallicity, and hence to opacity. To clarify the results obtained 
by Lenz et al. (2008) we compared OPAL and OP stellar models with the same 
mass (M=1.875 M Q , such as derived by Lenz et al. 2008) and at the same 
evolutionary stage (X c = 0.4). Their location in the HR diagrams is quite close, 
their radii differ by 0.05%, and hence their fundamental period, Ho, by 0.16% 
(log IIo=-0. 8380 instead of -0.8387). The effect on the period ratio is however 
larger going from 0.7697 for the OPAL model to 0.7722 for the OP one. This 
difference is an indication of a lower OP opacity with respect to OPAL values. 
In figure 1 we see that OP is lower than OPAL by a 10% in the region around 
logT - 6.05. 

As explained by Petersen & Jorgensen (1972), the period ratio (Ui/Uq) 
mainly depends on the effective politropic index (n = dlnP/dlnT — 1) of 
plasma in a region around x = r/R = 0.7, much deeper than the second ion- 
ization zone of helium. For their stellar envelope models this region corresponds 
to logT between 5.5 and 6.2. Figure 4, shows that the largest difference be- 
tween OP and OPAL opacity values is found close to x = 0.7. Different kr 
involves a different sound speed inside the star and hence different pulsation 
periods. A lower opacity implies larger values of period ratio. Hence, fitting 
with OP models the observational IIo and jTi/IIo of 44Tau required Lenz et 
al. (2008) to decrease by 10% the mass of the model, and therefore the ef- 
fective temperature and luminosity of the model were not compatible with the 
observational error box. 

Given the sensitivity of the period ratio to the opacity at logT <~ 6.05 
we wonder how the stellar parameters derived by using the Petersen's diagram 
depende on the uncertainty affecting the opacity values. We found that by 
increasing the OPAL opacities by at maximum a 5% around logT ~ 6.25, 
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a MS model with M = 1.9 M Q X = 0.70, Z = 0.02, X c = 0.195 and 
an overshooting parameter a = 0.3 satisfies the observational constrains (HR 
and Petersen diagrams locations) as shown in figure 5. This differs from the 
result found in Lenz et al. (2008), where only post-MS models were able to fit 
observations. 
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Figure 5: Left panel: evolutionary tracks for M = 1.9 M , Z = 0.02, X = 0.7, 
aov = 0.3 computed with OPAL opacity tables (solid line) and with OPAL enhanced 
by 5% in the region around logT ~ 6.05 (dashed line). The square corresponds to the 
observational error box (Lenz et al. 2008), and the star represent the model fitting 
IIo and the period ratio (IIi/IIo). Right panel: Petersen diagram corresponding to 
the evolutionary tracks in left panel. 



Other pulsators 
7 Doradus pulsators 

7 Doradus variables are F-type stars oscillating with high-order g-modes (Kaye 
et al. 1999), and the excitation mechanism has been identified by Guzik et 
al. (2000) as being due to the convective blocking of radiative flux. The loca- 
tion of the convective zone (and hence potentially stellar opacity) plays then a 
fundamental role in these pulsators. Computations considering the interaction 
pulsation-convection by Dupret et al. (2005) confirm the first results by Guzik 
obtained in the frozen convection approximation and shows that there is also 
an small contribution of the K-mechanism coming from the Z-bump of opacity 
at logT ~ 5.3, close to the bottom of the convective envelope in 7 Dor stars. 
Despite all that, changing opacity tables and the solar mixture do not signifi- 
cantly affect neither the instability strip nor the oscillation spectrum (Miglio et 
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al. in preparation). Nevertheless, these objects are still not well known, and 
some surprises could rise in the future. 

The Sun 

Oxygen is the main contributor to kr at the bottom of the solar convective 
envelope (logT ~ 6.3), followed by iron and neon. The decrease by 30% of 
solar metallicity resulting from Asplund et al. (2005) analysis has hence led to 
drop by 20% the opacity at the bottom of convective zone and to ruin the good 
agreement between the standard solar model and helioseismology (see Basu & 
Antia 2008 and references therein). The updated OP opacities imply only a raise 
of opacity at the bottom of the convective zone at maximum of 5%, too small 
to solve the problem of the helioseismic model. A large variety of suggestions 
have been proposed to recover the lost opacity and hence the agreement with 
helioseismic results, however none up to now (Basu & Antia 2008 for review) 
has provided satisfactory results, and the problem of the solar model remains. 
We would like, nevertheless, to recall that the good agreement between the 
Solar Standard Model (Christensen-Dalsgaard et al. 1996) and the seismic one 
was also damaged (even if the effect was small compared with the actual solar 
problem) when the updated version of OPAL opacity tables (Iglesias & Rogers 
1996) led to a decrease of opacity of 5% at the bottom of solar convective 
region logT ~ 6.3. The so called S-model was in fact computed with the first 
delivery of OPAL opacity tables (Rogers & Iglesias 1992). 

Conclusions 

In the nineties, the OPAL opacity tables signified a revolution in stellar physics. 
Since then, the progress in observational techniques and numerical compu- 
tations have highlighted new discrepancies between theory and observations. 
In spite of recent improvements some problems persist, for instance: i) the 
presence of (3 Cep pulsators in the SMC; ii) the excitation of high and low fre- 
quencies in B-type pulsators such as v Eri (Dziembowski & Pamyatnykh 2008); 
Hi) discrepancy between solar model and helioseismology; iv) the remaining 
15% disagreement between evolution and pulsation masses of cepheids vari- 
ables (Keller 2008 and references therein). Some of these discrepancies could 
be reduced by increasing (yes, again...) in limited domains of tempera- 
ture, in particular at logT > 6 (see also Zdrawkov & Pamyatnykh 2008, these 
proceedings). 

In spite of recent updates, some differences between OP and OPAL remains, 
mainly due to the different equation of state used in those approaches. The 
differences may seem small (~ 10 %), but we showed here that an enhancement 
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of 5% at logT ~ 6 may be able to modify the results of a seismic analysis. 
Both opacity compilations, OP and OPAL, were computed assuming reasonable 
approximations for the then current accuracy of observations. So, perhaps 
there is some room to increase stellar opacity. In fact, not all the elements are 
considered with the same degree of accuracy, and heavy and low abundance 
elements have not been taken in consideration since their contributions to kr 
was estimated to be small (Iglesias et al. 1995) for the solar model, and the 
computation very time consuming. 

Stellar pulsation theory has come out to be a powerful tool to probe different 
aspects of the stellar physics. Today it suggest that perhaps it is time to come 
back to review the basic input physics in stellar modeling, in particular stellar 
opacity computations. 

Acknowledgments. JM. and A.M. acknowledge financial support from the 
European Helio- and Asteroseismology Network HELAS, from the Prodex-ESA 
Contract Prodex 8 COROT (C90199) and from FNRS. 
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